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ABSTRACT 



Aims. We investigate the magnetic field power spectrum in the cool core galaxy cluster A2199 by analyzing the polarized emission of 
the central radio source 3C 338. 

Methods. The polarized radiation from the radio emitting plasma is modified by the Faraday rotation as it passes through the magneto- 
ionic intracluster medium. We use Very Large Array observations between 1665 and 8415 MHz to produce detailed Faraday rotation 
measure and fractional polarization images of the radio galaxy. We simulate Gaussian random three-dimensional magnetic field mod- 
els with different power-law power spectra and we assume that the field strength decreases radially with the thermal gas density as 
ril. By comparing the synthetic and the observed images with a Bayesian approach, we constrain the strength and structure of the 
magnetic field associated with the intracluster medium. 

Results. We find that the Faraday rotation toward 3C338 in A2199 is consistent with a magnetic field power law power spec- 
trum characterized by an index ^=(2.8+ 1.3) between a maximum and a minimum scale of fluctuation of A max =(35+28) kpc and 
A m in=(0.7+0.1) kpc, respectively. By including in the modeling X-ray cavities coincident with the radio galaxy lobes, we find a mag- 
netic field strength of (Z?o)=(11.7+9.0)yi/G at the cluster center. Further out, the field decreases with the radius following the gas 
density to the power of 77=(0.9+0.5). 
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1. Introduction 

In the past decade, significant progress has been made in char- 
acterizing the properties of the intracluster magnetic field. Most 
of what we know about intracluster magnetic field strength and 
structure is derived from the study of diffuse synchrotron emis- 
sion (radio halos, mini-halos, and relics) and Faraday rotation 
measures of polarized radio galaxies located inside or behind 
galaxy clusters (e.g., Carilli & Taylor 2002; Govoni & Feretti 
2004; Ferrari et al. 2008, and references therein). 

The study of the Faraday rotation of cluster radio galax- 
ies provides a detailed view of the intracluster magnetic field 
on scales < 100 kpc. Linearly polarized radiation propagating 
through a magnetized plasma experiences a rotation of the plane 
of polarization that is proportional to the thermal gas density 
and the magnetic field strength along the line-of-sight (e.g., Burn 
1966; Taylor et al. 1999). Indeed, it is possible to obtain impor- 
tant information about the intracluster magnetic fields by com- 
bining polarization images of radio sources located inside or be- 
hind galaxy clusters with X-ray observations of the thermal gas. 

The rotation measure (RM) distributions seen over extended 
radio sources are generally patchy, indicating that the intraclus- 
ter magnetic fields are not regularly ordered, but instead have 
turbulent structures on linear scales as small as a few kpc or 
less, both in merging and in relaxed clusters. Actually, magnetic 
field fluctuations seem to cover a wide range of spatial scales. 
In a few galaxy clusters containing radio sources with very de- 
tailed rotation measure images, the magnetic field power spec- 



trum has been estimated (EnBlin & Vogt 2003; Vogt & EnBlin 
2003; Murgia et al. 2004; Vogt & EnBlin 2005; Govoni et al. 
2006; Guidetti et al. 2008; Laing et al. 2008; Bonafede et al. 
2010; Guidetti et al. 2010; Kuchar & EnBlin 201 1). The rotation 
measure studies generally show central magnetic field strengths 
of a few fiG in merging galaxy clusters. Higher values of about 
10-40 yuG are typical of relaxed cooling core clusters (e.g. Dreher 
et al. 1987; Allen et al. 2001; Taylor et al. 2002), where the ob- 
served extreme RM magnitudes appear to be roughly propor- 
tional to the cooling flow rate (Taylor et al. 2002). In this paper, 
we investigate the magnetic field power spectrum in the nearby 
galaxy cluster A2199 by analyzing multifrequency polarization 
observations of the central radio source 3C338 taken with the 
Very Large Array (VLA). 

Chandra observations by Johnstone et al. (2002) revealed 
X-ray cavities in the cluster center associated with the radio 
lobes. The temperature of the intracluster medium decreases 
from 4.2 keV to 1 .6 keV over radii from 100 kpc to 5 kpc, imply- 
ing a drop in the radiative cooling time from 7 Gyr to 0.1 Gyr. 
These features seem to be consistent with a cooling flow, even 
if the action of radiative cooling should cool the central ICM 
to &r<lkeV. Therefore, as in many other cool core clusters, 
in A2199, it has been proposed that some heating mechanism 
(e.g., Tucker & Rosner 1983; Gaetz 1989; David et al. 2000) pre- 
vents the ICM from cooling down to temperatures below 1 keV. 
Kawano et al. (2003) suggest that previous AGN activity could 
be responsible for such heating in A2199. 
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Table 1. Details of the VLA observations of 3C 338. 



Obs. frequency, v 


Bandwidth 


VLA Configuration 


Time 


Data 


Project 


(MHz) 


(MHz) 




(h) 






1665 


50 


A 


8.5 


91-Jun-18/19 


GG0005 


4585/4685/4885 


50 


B 


7.0 


89-Apr-03 


AG0269 


4885 


50 


B 


1.5 


87-Dec-21 


AS0309 


4985 


50 


C 


3.0 


94-Nov-17 


BG0012 


4985 


50 


AB 


3.5 


95-Sep-ll/12 


BV0017 


8415 


50 


C 


3.0 


94-Nov-17 


BG0012 


8415 


50 


BC 


3.8 


00-Feb-26 


GG0038 



Col. 1: Observing frequency; Col. 2: Observing bandwidth; 

Col. 3: VLA configuration; Col. 4: Time on source; Col. 5: Dates of observations; Col. 6: VLA project name. 



Table 2. Relevant parameters of the total intensity and polarization images between 1665 and 8415 MHz. 



y 


Beam 


o-(I) 


cr(Q) 


cr(U) 


Sy 


(MHz) 


(" x " ) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy) 


1665 


2.5x2.5 


0.25 


0.029 


0.029 


1490+ 40 


4585 


2.5x2.5 


0.06 


0.033 


0.032 


520+ 20 


4685 


2.5x2.5 


0.06 


0.025 


0.026 


470+ 10 


4885 


2.5x2.5 


0.06 


0.026 


0.027 


440+ 10 


4985 


2.5x2.5 


0.08 


0.037 


0.037 


400+ 10 


8415 


2.5x2.5 


0.02 


0.015 


0.015 


170+ 10 



Col. 1: Observation frequency; Col. 2: FWHM; Col. 3,4,5: RMS noise of the I, Q, U images; 
Col. 6: Flux density. 



The radio source 3C 338 (otherwise known as B2 1626+39) 
is associated with the multiple nuclei cD galaxy NGC 6166, the 
brightest galaxy at the center of the cluster (Burns et al. 1983, 
Fanti et al. 1986), and it is classified as a restarting Fanaroff- 
Riley type I radio galaxy (see Murgia et al. 201 1 and references 
therein). On parsec scales, 3C 338 was the first established exam- 
ple of a symmetric, two-sided source in a radio galaxy (Feretti 
et al. 1993), consisting of a compact core and two symmetric 
relativistic jets (Giovannini et al. 2001). On kiloparsec scales, 
Burns et al. (1983) and Giovannini et al. (1998) showed that the 
radio structure consists of an active region, which includes the 
core and two symmetric jets terminating in two faint hot spots 
and two steep- spectrum radio lobes. The radio lobes are clearly 
associated with cavities in the central X-ray emission, as shown 
by Johnstone et al. (2002) and Gentile et al. (2007). Burns et al. 
(1983) pointed out a displacement between the large-scale struc- 
ture with respect to the restarted symmetric jets, which could 
indicate a possible motion of the central AGN inside the galaxy, 
but we speculate that the displacement could also be due to an 
interaction between the old radio lobes with bulk motions in the 
surrounding medium caused by the sloshing of the cluster core 
(Markevitch & Vikhlinin 2007). 

A2199 is an interesting target for Faraday rotation studies 
because the presence of X-ray cavities associated with the radio 
galaxy lobes indicates that the rotation of the polarization plane 
is likely to occur entirely in the intracluster medium, since com- 
paratively little thermal gas should be present inside the radio- 
emitting plasma. A previous rotation measure study has been 
done by Ge & Owen (1994) based on 5000 MHz VLA data. 
These authors show that 3C338 radio emission suffers a sig- 
nificant depolarization and is characterized by high RM values. 
By combining the information from this rotation measure im- 
age with deprojected ROSAT data, and assuming a very simple 
magnetic field model, Eilek & Owen (2002) infer an averaged 
magnetic field value along the line of sight of 15 yuG. 

In this paper we try to improve upon the previous estimate by 
analyzing additional data and by performing a numerical model- 



ing of the intracluster magnetic field fluctuations. We use VLA 
observations between 1665 and 8415 MHz to produce detailed 
Faraday rotation measure and fractional polarization images of 
the radio galaxy. Following Murgia et al. (2004), we simulate 
Gaussian random three-dimensional magnetic field models with 
different power-law power spectra, and we compare the synthetic 
and the observed images in order to constrain the strength and 
structure of the magnetic field associated with the intracluster 
medium. 

In § 2 we present the radio observations and the data reduc- 
tion. In §3 we describe the polarization properties of 3C338. 
In § 4 we present the Faraday rotation modeling. In sections 5 
and 6 we describe the results of the two-dimensional and three- 
dimensional simulations. Finally, in § 7 we summarize our con- 
clusions. Throughout this paper we adopt a ACDM cosmology 
with H = 71km s^Mpc" 1 , D m = 0.27, and Q A = 0.73. At 
the distance of 3C 338 (z=0.0311, Smith et al. 1997), 1" corre- 
sponds to 0.61 kpc. 
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Fig.l. Total intensity radio contours of 3C338 at 8415 MHz 
with an FWHM of 2.5" x 2.5". The first contour level is drawn at 
0.06 mJy/beam and the rest are spaced by a factor V2. The sensi- 
tivity (lcr) is 0.02 mJy/beam. The contours of the radio intensity 
are overlaid on the DSS2 red plate. 
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Table 3. Statistics of the rotation measure distribution. 



3 




Component 


(RM) 


CRM 


1 RM max | 


Err m 


(rm z ) i/z 




(rad/m 2 ) 


(rad/m 2 ) 


(rad/m 2 ) 


(rad/m 2 ) 


(rad/m 2 ) 


Total 


-54 


460 


903 


18 


463 


E-Lobe 


66 


473 


903 


20 


478 


W-Lobe 


-251 


356 


581 


16 


436 



RIGHT ASCENSION (J2000) 



Col. 1: Component of the source; Col. 2: Mean value of the RM distribution; 

Col. 3: Standard deviation of the RM distribution; 

Col. 4: Maximum absolute value of the RM distribution; 

Col. 5: Mean value of the RM fit error; 

Col. 6: Total power, (RM 2 ) 1/2 = ((RM) 2 + cr 2 ^) 1 ! 2 . 
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3C 338 8415 MHz 
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Fig. 2. Total intensity contours and polarization vectors at 1665 
(top), 4585 (middle), and 8415 MHz (bottom). The angular res- 
olution is 2.5" x 2.5". Contour levels start at 3<xi, and increase 
by factors of two. The lines give the orientation of the electric 
vector position angle (E-field) and are proportional in length to 
the fractional polarization (1" ~ 5%). 



figuration and frequency, all visibilities corresponding to the 
same frequency were combined using the task DBCON in order 
to improve (u, v) coverage and sensitivity. The combined data 
were then self-calibrated to produce the final images. The im- 
ages of total intensity / and Stokes parameters Q and U were 
restored with a circular beam of 2.5". The noise levels of /, 
U, and Q are summarized in Table [2l Polarized intensity (cor- 
rected for the positive bias), P = ^Q 2 + U 2 , polarization angle, 
Y = 0.5tan~ l (U/Q), and fractional polarization FPOL = P/I 
images were produced for each frequency. 

The total intensity contours of the radio galaxy 3C 338 over- 
laid on the DSS2 red plateQ are showed in Fig.HJ The 8415 MHz 
radio galaxy emission peaks at RA(J2000)=16 h 28 m 38 s and 
DEC(J2000)=+39 d 33 m 04 s . It is also the same location of the X- 
ray peak to within the image resolution. Therefore we chose this 
reference location for our analysis. 

3. Polarization analysis 

In Fig. [2] we present the total intensity radio contours and po- 
larization vectors at 1665, 4585, and 8415 MHz at a resolution 
of 2.5" x 2.5". The length of the vectors is proportional to the 
fractional polarization, while their orientation is the same as the 
projected E-field. The fractional polarization and the polariza- 
tion angle were obtained by considering only those pixels where 
the fractional polarization is above 3<xfpol at the corresponding 
frequency. The total intensity emission of the source extends out 
to about V from the cluster center. The total flux density of the 
radio source and all the relevant parameters of the /, g, and U 
images are reported in Tabled 



2. Radio observations and data reduction 

We present polarimetric archival observations of 3C 338 at 1665, 
4585, 4685, 4885, 4985, and 8415 MHz performed at the VLA 
in the A, B, and C configurations between 1987 December and 
2000 February. All observations were made with a bandwidth 
of 50 MHz. The pointing center is identified by the coordinates 
RA(J2000)=16 h 28 m 38 s , and DEC(J2000)=+39 d 33 m 04 s , and the 
details are provided in Table[T] 

All data were reduced using the NRAO's Astronomical 
Image Processing System (AIPS) package. The radio source 
3C286 was used both as primary flux density calibrator and 
as reference for the absolute polarization angles. Phase cali- 
bration was derived from nearby sources, periodically observed 
over a wide range in parallactic angle to separate the source 
polarization properties from the antenna polarizations. We pro- 
ceeded with the standard calibration procedures, imaging and 
self-calibration. Once all the data were calibrated for each con- 



3.1. Rotation measure 

The presence of a magnetic field in an ionized plasma creates a 
difference in the phase velocities for left versus right circularly 
polarized radiation. As a consequence, the polarized emission 
from a radio source propagating through the plasma experiences 
a phase shift between the two components. This corresponds to 
a rotation in the polarization angle. For a completely foreground 
screen (Burn 1966), which we expect here, the rotation is 

¥obsW = ¥int + ^RM, (1) 

where ^obsW) is the observed position angle at wavelength A, 
^int the intrinsic polarization angle of the polarized emission, 
and RM the rotation measure. Under this assumption, by observ- 
ing radio galaxies at different wavelengths, the rotation measure 
can be derived from a linear fit of the polarization position angle 
^obsW) versus A 2 . By considering an electron density n Q (cm -3 ), 

1 http://archive.eso.org/dss/dss 
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a magnetic field B (juG), and a path length 1 (kpc), the Faraday 
RM is 

RM = 812 B^g] -dl[k P c] radm 2 . (2) 

Jo 

We produced the Faraday RM image by running the 
FARADAY code (Murgia et al. 2004). The software requires as 
input Q and U images for each frequency and outputs the RM, 
the intrinsic polarization angle their errors images, and^ 2 
maps. The RM image is created pixel by pixel by fitting the ob- 
served polarization angle ^obs versus the squared wavelength A 2 
for all the frequencies. To reduce the problems associated with 
nn ambiguities, the fitting algorithm can perform a sequence of 
improvement iterations. In the first iteration, only a subset of 
high signal-to-noise pixels is considered. In the successive itera- 
tions, lower signal-to-noise pixels are gradually included and the 
information from the previous iteration is used to assist the fit of 
the A 2 law. 

We considered only those regions of the radio source where 
the total intensity emission at 8415 MHz is above 5<xt. Only pix- 
els with an uncertainty in polarization angle below 10° at each 
wavelength were considered. The RM fit is generated if this con- 
dition is satisfied for at least five frequency maps. Almost half of 
the total number of the pixels in the RM image are based on a 
five-frequency fit. The remainder are based on a six-frequency 
fit. 

The final rotation measure image of 3C338 is shown in 
Fig. [3] (top left panel) with total intensity contours at 8415 MHz 
overlaid. The image has a resolution of 2.5", which corresponds 
to 1.5 kpc at the distance of 3C 338. The RM has a patchy struc- 
ture with values ranging from -1300 to 900rad/m 2 . As shown 
by the histogram in Fig. [3] (top right panel), the RM distribution 
is characterized by a mean value of (RM) = -54rad/m 2 
and a standard deviation cr RM = 460rad/m 2 . The RM fit is 
characterized by a mean error Err m = 1 8 rad/m 2 . In Tabled we 
report the statistics of the RM distribution of the entire source 
and, separately, of the east and west lobes. The mean and 
dispersion of the RM are quite different, but the "total power" 
{RM 2 ) 112 = ((RM) 2 + cTr M ) 1/2 is very similar. Therefore, we 
cannot constrain the inclination to the line of sight of the source, 
and in the following analysis we assume that 3C338 is on the 
plane of the sky. Examples of the observed position angle ^obs 
as a function of A 2 for four high signal-to-noise pixels are shown 
in the bottom panels of Fig.[3j The observed fits are linear, as 
expected in the case of a foreground Faraday screen. 

We have estimated the contribution of our own Galaxy to the 
Faraday rotation in the direction of 3C 338. In galactic coordi- 
nates, the radio source is located at / = 62.9° and b = 43.7°. 
We computed the average of the RM values reported in Taylor 
et al. (2009) for a region of about 3° around this direction and 
from this estimated a Galactic contribution of -13 rad/m 2 . Since 
the Galactic foreground is small compared to the RM intrinsic to 
3C 338, it has a negligible impact on the RM. 

3.2. Depolarization 

In Table|4]we report the average fractional polarization levels of 
3C338 obtained by considering the same pixels as we used to 
calculate the RM image. The radio source is less polarized at 
longer wavelengths. This behavior can be interpreted in terms of 
variations is the RM on smaller scales than the beam of the radio 
images. 



Table 4. Fractional polarization between 1665 and 8415 MHz 



y 


FPOL 


(MHz) 


(%) 


1665 


1.1+0.3 


4585 


13.6+0.3 


4685 


15.8+0.2 


4885 


18.2+0.3 


4985 


23.0+0.4 


8415 


41.7+0.6 



Col. 1: Observation frequency; 
Col. 2: Fractional polarization. 



These unresolved RM structures in the foreground screen 
cause a depolarization of the signal which, to first order, can be 
approximated by 

FPOL = FPOL exp(-aA 4 ), (3) 

where a is a value related to the RM gradient within the ob- 
serving beam and FPOLq is the intrinsic fractional polarization 
(Burn 1966, see also Laing et al. 2008 for a more recent deriva- 
tion). We fitted the Burn law only between 4585 and 8415 MHz 
(6 and 3.6 cm) observations, since the formula is not applicable 
in the long-wavelength regime (20 cm/1665 MHz) (see Tribble 
1991 and g5]). We find a=(66+6)xl0 3 rad 2 /m 4 . 



4. Faraday rotation modeling 

Our aim is to constrain the intracluster magnetic field strength 
and structure in A2199 by using the information from the ra- 
dio galaxy RM and polarization images presented in the previ- 
ous section. Our modeling is based on the assumption that the 
Faraday rotation is occurring entirely in the intracluster medium. 
In particular, we suppose that there is no internal Faraday rota- 
tion inside the radio lobes, as suggested by the X-ray cavities 
detected by Johnstone et al. (2002) and by the observed linearity 
of the polarization angle ^obs versus A 2 (see Fig. [3]). Moreover, 
we suppose that any possible local RM enhancement occurring 
at the interface between the radio lobes and the surrounding 
medium is negligible compared to the total RM across the cluster 
(but see Rudnick & Blundell 2003 for a contrary viewpoint). 

As shown in Eq.O the rotation measure is the integral along 
the line-of- sight of the intracluster magnetic field and the ther- 
mal gas density. To derive useful information about the intraclus- 
ter magnetic field we first need a model for the spatial distribu- 
tion of the thermal electron density. 

4.1. Thermal gas modeling 

The distribution of the thermal electrons was modeled using the 
Chandra X-ray observation of A2199 by Johnstone et al. (2002). 
In the left panel of Fig. [4] we show a superposition of the VLA 
8415 MHz contour levels on the Chandra X-ray image in the 
0.2-10 keV band of the central region of A2199. In the top right- 
hand panel of Fig. [4] we show the deprojected thermal gas den- 
sity n Q (r) profile derived by Johnstone et al. (2002). In the bottom 
righthand panel we show the rotation measure pixel distribution 
as a function of the distance from the cluster center. The RM 
pixels are sampled at distances from the cluster core in the range 
3 to 30 kpc, with a higher concentration between 10 and 20 kpc. 
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Fig. 3. Top left: total intensity radio contours of 3C 338 at 8415 MHz overlaid on the RM image. The angular resolution is 2.5" x2. 5". 
Contour levels are drawn at: 0.06, 0.12, 0.24, and 0.96mJy/beam. Top right: the histogram of the RM values. Bottom: sample plots 
of the observed position angle versus the squared wavelength at different source locations. The inset is a zoom of the observed 
polarization angle corresponding to frequencies between 4585 and 8415 MHz. 
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1 2 3 4 5 6 7 




RIGHT ASCENSION (J2000) r [kpc] 

Fig. 4. Left panel: Adaptively smoothed Chandra X-ray image (Obs. ID 497) in the 0.2-10 keV band of the galaxy cluster A2199 
with the VLA contour levels at 8415 MHz overlaid. Pronounced X-ray cavities can be identified in correspondence of the radio 
galaxy lobes. Top right panel, deprojected thermal gas density profile (Johnstone et al. 2002). Bottom right panel, pixel distribution 
of the rotation measure image as a function of the distance from the cluster center. 



We initially fit the deprojected thermal gas density profile 
with a single /3-model 



n Q (r) = n 1 + — 



(4) 



where r is the distance from the cluster X-ray centroid, no the 
central electron density, and r c the cluster core radius (Cavaliere 
& Fusco-Femiano 1976). But this thermal gas model does not al- 
low a good description of the central jump owing to the high cen- 
tral density, typical of cooling core galaxy clusters. Therefore, 
we tried to improve the fit by considering a double /3-model, 



n Q (r) = n 0in 



r 2 



1 + — 

ri 



+ n 0sxt 1 + -r- 



r 2 
rl 



(5) 



where no int and no ext are the internal and external central electron 
densities, while r Cint and r Cext are the internal and external cluster 
core radii, respectively. We do not discuss the physical validity 
of this model here. Our purpose is to empirically obtain an ana- 
lytic profile that provides a reasonable continuous description of 
the deprojected gas density observed by Johnstone et al. (2002). 
This analytic profile can be conveniently extrapolated to both the 
cluster center and to large radii so that the integral in Eq.[2l can 
be performed. 

The double /3-model best-fit parameters are reported in 
TableO while the best-fit profile is shown as a continuous line 
in the righthand panel of Fig.|4l In the same plot we also show, 
as a dotted line, the external part of the model that has been con- 
strained by a fit of the data points with r > 18 kpc. The total 
central density we find is no = fto int + fto ext =0.101 cm -3 . 



Table 5. Double J3- model parameters. 



PARAMETER 


VALUE 


UNITS 


Ant 


1 5 +o:1 

-i>-5 




r q n t 


9+2 


kpc 




074+ - 004 


cm -3 


A* 


39+ 01 

-01)3 




r Cext 


26 - 0+ 60 


kpc 


^0 ext 


027 +o;o03 

7 -0.003 


cm -3 



4.2. Magnetic field modeling 

The patchy structure characterizing the RM image in Fig. [3] can 
be interpreted as a signature of the turbulent intracluster mag- 
netic field. In particular, the dispersion of the RM values can be 
related to the strength and structure of the magnetic field. If the 
magneto-ionic medium is approximated by uniform cells of size 
A c with random orientation in space, the Faraday rotation from 
a physical depth L » A c is expected to be a Gaussian with zero 
mean and dispersion given by 



a 2 = (RM 2 ) = 812 2 A ( 



f (n { 
Jo 



Q B z ) 2 dl, 



(6) 



where B z is the magnetic field component along the line-of- sight 
(e.g., Felten 1996). 

More generally, the magnetic field auto-correlation function 
Ce z (r) and the rotation measure auto-correlation function Crm 
are related by the Abel transform (EnBlin & Vogt 2003), 



RM 



(r ± ) = 812 2 ^Lx 



■1 



dr 



(7) 
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Parameter 


Description 


Means of investigation 


(Bo) 


strength at the cluster center 


3D simulations 


V 


radial slope, (B(r)) = (B ) (*®f 
power spectrum index, |Z? k | 2 oc k~ n 
minimum scale of fluctuation, A min = 27r/£ max 


3D simulations 


n 

A min 


2D simulations 
2D simulations 


A ma x 


maximum scale of fluctuation, A max = 2^/^ min 


2D simulations 



where, for simplicity, the thermal gas density n Q is considered 
constant here. 

If the characteristic fluctuation scale of the magnetic field is 
much smaller than the physical depth of the Faraday screen (such 
that (B) = and (RM) = 0), in the limit r ± -> 0, we obtain 



°~RM 

where 



812 2 ^LA B oi z 



2jT L/2 C Bz (r>/r 
C Bz (0) 



(8) 



(9) 



defines the magnetic field auto-correlation length, which should 
be the appropriate scale to use for A c in Eq.|6] 

For an isotropic divergence-free random field, it can be 
shown that A B is related to the magnetic field power spectrum, 
|5kl 2 , by 



3n £*\B k \ 2 kdk 
2 JT \B k \Wdk 



(10) 



where k = 2n/A is the wave number (EnBlin & Vogt 2003). 

This is a particularly important result, since it implies that 
we must know the power spectrum of the magnetic field fluctu- 
ations to determine the strength of the field from the RM image. 
We also note that the RM power spectrum is proportional to the 
magnetic field power spectrum, 



\RM k \ 2 oc nlL\B k \\ 
but the RM auto-correlation length, 
J °° \RM k \ 2 dk 



Ai 



RM 



' f °° \RM k \ 2 kdk 



(11) 



(12) 



is not the same as A B . 

To limit the number of free parameters, in this work we 
choose to model a power-law power spectrum^ with index n of 
the form 



\B k \ 2 oc r n 



(13) 



in the wave number range from k m { n to & max and outside. 
Moreover, we suppose that the power- spectrum normalization 
varies with the distance from the cluster center such that the av- 
erage magnetic field strength scales as a function of the thermal 
gas density according to 



(B(r)) = (B ) 



n 



(14) 



2 Throughout this paper the power spectra are expressed as vectorial 
forms in &-space. The one-dimensional forms can be obtained by mul- 
tiplying by Ank 2 and Ink the three and two-dimensional power spectra, 
respectively 



where (Bo) is the average magnetic field strength at the center 
of the cluster, and n Q (r) is the thermal electron gas density. We 
have to note that the tapering of the magnetic field power spec- 
trum due to the thermal gas density may not preserve exactly 
the magnetic-field power- spectrum power-law shape at the edges 
(see Fig.[TO|. 

Overall, our magnetic field model depends on the five pa- 
rameters listed in Tabled the strength at the cluster center (Bo), 
the radial slope rj, the power spectrum index n, and finally the 
minimum and maximum scales of fluctuation, A m i n = 2n/k max 
and A max = 2n/k m i n , respectively. 

4.3. Bayesian inference 

To constrain the magnetic field strength and structure, we pro- 
ceeded in two steps. First, we performed a two-dimensional anal- 
ysis of the RM fluctuations and of the source depolarization 
to constrain the slope n and the range of scales of the power 
spectrum. Second, we performed three-dimensional numerical 
simulations to constrain the strength of the field and its scal- 
ing with the gas density. In both cases, we made use of the 
FARADAY code to produce synthetic polarization images of 
3C338 and to compare them to the observed ones. We com- 
pared model and data using the Bayesian inference, whose use 
was first introduced in the RM analysis by EnBlin & Vogt (2003). 
Because of the random nature of the intracluster magnetic field, 
the RM image we observe is just one possible realization of the 
data. Different realizations of magnetic field characterized by the 
same power spectrum will generate different RM images. Thus, 
rather than try to determine the particular set of power spectrum 
parameters that best reproduces the given realization of the data, 
it is perhaps more meaningful to search for that distribution of 
model parameters that maximizes the probability of the model 
given the data. The Bayesian inference offers a natural theoreti- 
cal framework for this approach. 

The Bayesian rule relates our prior information on the distri- 
bution P(0) of model parameters to their posterior probability 
distribution P(6 | D) after the data D have been acquired 



P(0\D) 



UP | 0)P(0) 
P(D) ' 



(15) 



where L(D | 0) is the likelihood function, while P(D) is called 
the evidence. The evidence acts as a normalizing constant and 
represents the integral of the likelihood function weighted by 
the prior over all the parameters space 



P(D) 



I 



L(D | 0)P(0)d0. 



(16) 



The most probable configuration for the model parameters is 
obtained by maximizing the joint posterior given by the product 
of the likelihood function with the prior probability. We used a 
Markov Chain Monte Carlo (MCMC) method to extract sam- 
ples from the posterior probability distribution. In particular, we 
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implemented the Metropolis-Hastings algorithm, which is capa- 
ble of generating a sample of the posterior distribution without 
the need to calculate the evidence explicitly, which is often ex- 
tremely difficult to compute since it would require exploring the 
entire prior space. The MCMC is started from a random initial 
value #o and the algorithm is run for many iterations by selecting 
new states according to a transitional kernel, Q(0, 0), between 
the actual, 6, and the proposed position, 0' . The proposed posi- 
tion is accepted with probability 



h = min 



L(D\0')P(0')Q(0',0) 



L(D\0)P(0)Q(0,0>) 



(17) 



We chose for Q a multivariate Gaussian kernel. The MCMC 
starts with a number of "burn-in" steps during which, according 
to common practice, the standard deviation of the transitional 
kernel is adjusted so that the average acceptance rate stays in 
the range 25% -50%. After the burn-in period, the random walk 
forgets its initial state and the chain reaches an equilibrium dis- 
tribution. The burning steps are discarded, and the remaining set 
of accepted values of is a representative sample of the posterior 
distribution that can be used to compute the final statistics on the 
model parameters. 

5. 2D analysis 

We performed a preliminary two-dimensional analysis of the 
RM fluctuations and of the source depolarization to constrain 
the slope n and the range of scales of the power spectrum. 
The two-dimensional analysis relies on the proportionality be- 
tween the magnetic field and the rotation measure power spec- 
tra. On the basis of this proportionality, the index n of the two- 
dimensional rotation measure power spectrum is the same as the 
three-dimensional magnetic field power spectrum: 

|2 oc r n . 



\RM k 



(18) 



At this stage, we are only interested in the shape of the power 
spectrum and not in the exact value of the normalization. Indeed, 
as a first approximation, we do not consider the effect of the spa- 
tially variable gas density; i.e., we suppose that the depth of the 
Faraday screen is on average the same over all of the source. 
We simulated synthetic images with a given power spectrum 
in a two-dimensional grid. The simulations start in the Fourier 
space, where the amplitudes of the RM components are selected 
according to Eq.[l8l while the phases are completely random. 
The RM image in the real space is obtained by a fast Fourier 
transform (FFT). 

As a compromise between computational speed and spa- 
tial dynamical range, we adopted a computational grid of 
2048x2048 pixels with a pixel size of 0.05 kpc. This grid al- 
lowed us to explore RM fluctuations on scales as small as A m i n ^ 
0.1 kpc, i.e. one order of magnitude smaller than the beam of 
the observations (FWHM 1.5 kpc). Simultaneously, we were 
able to investigate fluctuations as large as A max ^ 100 kpc, i.e., 
comparable to the linear size of 3C 338. We determined the dis- 
tribution of the best-fit parameters using the Bayesian approach 
outlined in 34.31 In particular, we characterize the RM image by 
its structure function 

WrJ = ([RM(x,y) - RM(x + Ax,y + Ay)]\ x , yh (19) 

which is obtained by averaging the difference in RM values 
corresponding to pixels located at the scale distance r ± = 
^jAx 2 + Ay 2 , and is related to the RM auto-correlation function 



by the simple relation 

«rj = 2(cr 2 M + (RM) 2 ) - 2C RM (rJ. 



(21) 



The observed structure function is shown in the top lefthand 
panel of Fig.0 the formal error bars are comparable to the size 
of the dots. For an isotropic field, the RM auto-correlation func- 
tion and power spectrum are related by the Hankel transform 



Crm(t±) 
where 



" 2 "l 



J (krJ\RM k \ 2 kdk, 



Jo(kr ± ) 



2?r Jo 



-ikr ± cos# 



de 



(22) 



(23) 



is the zero-order Bessel function. For the power-law power 
spectrum in Eq.QjO we can identify three regimes in the RM 
structure function. In the asymptotic small- separation regime, 
r ± <?c 2n/k max , the structure function increases as Srm^x) 00 
In the intermediate regime, where 2n/k mSiX «: r ± <?c 2n/k m i n , 
SrmOx) oc r n ~ 2 . Finally for r ± » 2n/k m i n , the structure function 
saturates to the constant value of Srm - 2<x 2 , M . However, it can 
be hard to discern these three regimes directly from the observed 
structure function because of the coarse resolution of the radio 
images and the undersampling of the large separations because 
of the finite size of the RM image. Indeed, we need to resort to 
numerical simulations in order to account for the effects of these 
window functions. The synthetic images are thus gridded to the 
same geometry as the data and are convolved to the same angu- 
lar resolution. Moreover, we masked the synthetic images using 
the observations to reproduce the window function imposed by 
the shape of 3C 338, and we added Gaussian noise with an rms 
value of Err$x, in order to mimic the noise of the observed RM 
image. 



Crm(^) = (RM(x,y) • RM(x + Ax, y + Ay)> (xo 



(20) 



We then applied the Bayesian method described in § !4.3l bv 
choosing uniform priors for the four power spectrum parame- 
ters: the normalization norm, the slope n, the minimum, and the 
maximum wave number k m m and k max . The statistics of the sim- 
ulated RM fluctuations are Gaussian, so their structure function 
at a given separation has a log-normal distribution. At each step 
of the MCMC, we performed 30 realizations of the same power 
spectrum with different random seeds and evaluated the aver- 
age and the dispersion of the logarithm of the model structure 
function. We averaged the offsets and the dispersions at all the 
individual separations. From these values we computed the like- 
lihood of the observed structure function. 

In Fig. [5] we show the results of the Bayesian analysis of the 
rotation measure structure function. The output from the MCMC 
is a 4-dimensional hypercube containing a sample of the pos- 
terior distribution of the power spectrum parameters. The top- 
right, bottom-left, and bottom righthand panels show the two- 
dimensional marginalization of this hypercube as color images. 
In addition, one-dimensional marginalizations are shown as his- 
tograms along each axis of these images. The two-dimensional 
marginalizations represent the projected density of samples in 
the MCMC, which is proportional to the posterior probability of 
that particular couple of model parameters. The maximum scale 
of the fluctuation, the normalization, and the slope of the power 
spectrum appear to be characterized by a peak that corresponds 
to the maximum posterior probability for that configuration of 
parameters, while for the minimum scale of fluctuation we only 
have an upper limit. 

To provide a visual comparison between model and data, in 
the top lefthand panel of Fig. [5] we show the observed structure 
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Fig. 5. Top left panel. Bayesian analysis of the RM structure function. The dots represent the data (error bars are comparable to the 
size of the symbols). The shaded area represents the population of synthetic RM structure functions from the posterior distribution. 
The dashed line corresponds to the most probable value for the model parameters (see text). Top right, bottom left, and bottom right 
panel, one-dimensional (histograms) and two-dimensional (colors and contours) marginalization of the posterior for the model 
parameters. The contours are traced at 0.9, 0.75, and 0.5 of the peak value. 



function along with the population of synthetic structure func- 
tions contained in the posterior sample and the structure function 
for which posterior is maximum (best fit). The brighter pixels in 
the shaded image occur where many different synthetic struc- 
ture functions overlap each other. The probability of the model 
given the data is higher in the brighter regions and is lower 
in the darker regions. The shaded region is relatively narrow 
for small separations but widens significantly for r± > lOkpc, 
indicating that our sensitivity to the large-scale separations is 
lower. Overall, the data stay close to the high-probability re- 



gion, indicating that the model is doing a good job reproduc- 
ing the observed RM structure. However, for r± > 20kpc the 
observed structure function decreases while the model stays rel- 
atively constant. The turnover of the observed structure function 
is very likely due to the lower gas density at large distances from 
the cluster center, which results in a systematic decrease in the 
power spectrum normalization, hence of the Faraday rotation on 
large scales. This effect is not included in the current modeling, 
but it will be investigated with the aid of three-dimensional sim- 
ulations in §0 
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Fig. 6. Top left panel. Bayesian analysis of the source depolarization. The dots represent the data. The shaded area represents the 
population of synthetic polarization from the posterior distribution. The dashed line corresponds to the most probable value for 
the model parameters (see text). Top right, bottom left, and bottom right panel: one-dimensional (histograms) and two-dimensional 
(colors and contours) marginalization of the posterior for the model parameters. The contours are traced at 0.9, 0.75, and 0.5 of the 
peak value. 



The power spectrum used to model the RM image should 
also be consistent with the observed depolarization of the radio 
source at increasing wavelengths (see 93.21) . The depolarization 
is caused by RM variations on smaller scales than the beam that 
results in an incoherent sum of the radio signal. Indeed, mod- 
eling of the polarization amplitude can be used to place more 
stringent constraints on the minimum scale of fluctuation of the 
magnetic field and on the slope of the power spectrum. Using 
the same simulation set-up presented above, we reproduced the 
expected polarized signal from 3C 338 as a function of the wave- 



length. We first constructed an image of the source polarization 
at A = 0, which is characterized by an intrinsic degree of polar- 
ization FPOLq. We then simulated different RM images at full 
resolution and used them to rotate the intrinsic polarization vec- 
tors according to Eq.[T] These full-resolution images are finally 
convolved at the same resolution as the observations, resulting 
in beam depolarization of the signal at longer wavelengths. 

We used the Bayesian inference and a Gaussian likelihood 
to estimate the distribution of the power spectrum parameters, 
which maximizes the probability that the observed depolariza- 






tion is a realization of the model. In Fig.[6]the posterior from the 
depolarization analysis is shown for the five free parameters: the 
intrinsic degree of polarization FPOLq, the normalization norm, 
the slope n, the minimum, and the maximum scale of fluctua- 
tion A m i n , and A max of the magnetic field power spectrum. In 
the top right and bottom panels the two-dimensional (colors and 
contours) and one-dimensional (histograms) marginalization of 
the posterior are shown for three different combinations of the 
model parameters. All the model parameters appear to be well 
constrained, and their values are consistent with the structure 
function analysis. In the top lefthand panel, the observed frac- 
tional polarization as a function of the fourth power of the wave- 
length is shown, together with the sample of realizations from 
the posterior and the best fit. In Fig. [7] the observed fractional po- 
larization images are compared with synthetic realizations cor- 
responding to the best fit parameters. The synthetic fractional 
polarization trend at higher frequencies (4585 to 8415 MHz) is 
consistent with a Burn law, FPOL = FPOLq exp(-<2/i 4 ), with 
^synth=(61±3)xl0 3 rad 2 /rn 4 , in very good agreement with the ob- 



served value <2 o bs=(66±6)xl0 3 rad 2 /rn 4 . However, it is clear from 
the top lefthand panel of Fig.[6|that our simulations are also able 
to explain the observed polarization levels reasonably well at low 
frequencies (1400 MHz) where the Burn law breaks down. 

Overall, the combined two-dimensional analysis, RM image, 
and depolarization allowed us to obtain a first insight into the 
shape of the magnetic field power spectrum. We constrained 
the power spectrum index to ^=(2.8+1.3), and the minimum 
and maximum scales in the range from A m i n =(0.7±0.1)kpc to 
A ma x=(35±28)kpc, where the given errors represent the disper- 
sion of the one-dimensional marginalizations. In the next step 
we fix these values and constrain the strength of the magnetic 
field and its scaling with the gas density with the aid of three- 
dimensional simulations. 

6. 3D simulations 

We can construct a three-dimensional model of the intracluster 
magnetic field by following the numerical approach described 
in Murgia et al. (2004). The simulations begin in Fourier space 
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by extracting the amplitude of the magnetic field potential vec- 
tor, A(k), from a Rayleigh distribution whose standard deviation 
varies with the wave number according to |Ak| 2 oc k~ n ~ 2 . The 
phase of the potential vector fluctuations is taken to be com- 
pletely random. The magnetic field is formed in Fourier space 
via the cross product B(k) = ik x A(k). This ensures that the 
magnetic field is effectively divergence free. We then perform a 
three-dimensional FFT inversion to produce the magnetic field in 
the real space domain. The field is then Gaussian and isotropic, 
in the sense that there is no privileged direction in space for the 
magnetic field fluctuations. The power- spectrum normalization 
is set such that the average magnetic field strength scales as a 
function of the thermal gas density according to Eq.[l4l The gas 
density profile is described by the double-beta model in Eq.[5j 
Actually, we know that the gas density distribution in A2199 
deviates from spherical symmetry because of the X-ray cavities 
corresponding to the radio lobes of 3C 338. Hence, we include 
the X-ray cavities in the three-dimensional simulations by re- 
moving from the gas density modeling all the magneto-ionic ma- 
terial inside two ellipsoidal regions centered on the radio lobes, 
as shown in the top panel of Fig.[8j In the bottom panel we trace 
the profile of the X-ray brightness along an horizontal slice pass- 
ing through the radio lobes for the double beta-model with and 
without cavities. The model with cavities provides a better de- 
scription of the observed X-ray brightness along the slice. 

We simulated the random magnetic field by using a cubical 
grid of 1024 3 pixels with a cell size of 0.16kpc/pixel. The syn- 
thetic RM images were obtained by integrating numerically the 
gas density and the magnetic field product along the line-of- sight 
accordingly to Eq.[2 In a similar way to the two-dimensional 
simulations, the synthetic RM images were gridded, convolved, 
blanked, and noised as the observed RM image before the com- 
parison with the data. We used the Bayesian approach to find the 
posterior distribution for (Bo) and n, which maximizes the prob- 
ability that the observed structure function is a realization of the 
model. The slope n and the scales k m { n and & max were kept fixed 
at the values found with the two-dimensional analysis. 

We started with uniform priors for (Bo) and n and we 
evaluated the likelihood of the structure function by mean of 
30 different configurations for the magnetic field phases at 
each step of the MCMC. The result of the Bayesian analy- 
sis is shown in Fig. [9) In the bottom panel we present the 
two-dimensional (colors and contours) and one-dimensional 
marginalizations of the posterior. The two parameters appear 
well constrained. We found a magnetic field with a central 
strength <fl >=(11.7±9.0)//G, and a radial slope ?7=(0.9±0.5). 
In the top panel we show the observed structure function along 
with the synthetic structure functions from the posterior and the 
best fit. The three-dimensional modeling represents a signifi- 
cant improvement over the two-dimensional analysis. In fact, 
we are now able to describe the overall shape of the observed 
structure function with good accuracy, including the turnover at 
large separations that is most likely due to the decrease in the 
Faraday rotation with radius. In the three-dimensional simula- 
tions we kept the maximum scale of the magnetic field fluctu- 
ation fixed to A max =35 kpc, which is the value found from the 
two-dimensional analysis in §[5j The consequences of a differ- 
ent choice of this parameter are discussed in further detail in the 
Appendix where we also tested the cases A max =10 and 164 kpc. 
The value A max =35 kpc still provides the best description of the 
observed RM structure function. In Fig.[K)|we plot the magnetic 
field power spectra and radial profiles corresponding to the max- 
imum posterior for the three values of A max . 
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Fig. 8. Top panel: X-ray cavity model. The darker region indi- 
cates the location where the thermal gas has been subtracted 
from the double /3-model. Bottom panel: horizontal slice passing 
through the X-ray cavities at the cluster center. Continuous line 
represents the observations from the Chandra image, dashed line 
simulations with cavities, dotted line simulations without cavi- 
ties. 



We performed the same Bayesian analysis also not taking 
of the X-ray cavities into account (not shown). We found very 
similar results to those including the cavities and we conclude 
that the presence of these voids in the gas density distribution 
has a second-order impact on the Faraday rotation measures we 
are analyzing. This could be because the most of the lines of 
sight sampled from our observed RM image do not intercept the 
regions of the cluster affected by the cavities. 

Finally, in Fig.Qj] we present a qualitative comparison be- 
tween the observed and a synthetic RM images taken from 
the simulations with cavities. Although very simple, the power- 
spectrum model adopted here provides a reasonable description 
of the patchy RM structure seen in the data. 

It is interesting to compare our result with the independent 
estimate based on the theoretical work by Kunz et al. (201 1). By 
assuming that turbulent dissipation balances radiative cooling at 
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all radii inside the cluster core, Kunz et al. (2011) find in the 
bremsstrahlung regime (that is T > 1 keV) 

where T is the temperature, while £ is expected to range between 
0.5 and 1 in a turbulent plasma. By considering a central den- 
sity cm -3 and a central temperature T^2 keV, according 
Johnstone et al. (2002), they found ^o-l l^" 1/2 yuG, very close to 
our estimate based on the Faraday rotation. 

7. Conclusions 

We investigated the magnetic field power spectrum in the cool 
core galaxy cluster A2199 by analyzing the polarized emis- 
sion of the central radio source 3C 338. We used archival VLA 
observations between 1665 and 8415 MHz to produce detailed 
Faraday rotation measure and fractional polarization images of 
the radio galaxy. We observed a significant depolarization of the 
radio emission of the radio galaxy and high RM values. This 
agrees with the analysis performed by Ge& Owen (1994) at 
5000 MHz, although our result takes the additional information 
at 1665 and 8415 MHz into account. 

We simulated Gaussian random three-dimensional magnetic 
field models with different power-law power spectra, and we 
assumed that the field strength decreases radially as a power 
of the thermal gas density as n^. By comparing the synthetic 
and the observed images with a Bayesian approach, we con- 
strained the strength and structure of the magnetic field associ- 
ated with the intracluster medium. We found that the Faraday 
rotation toward 3C338 in A2199 is consistent with a mag- 
netic field power-law power spectrum characterized by an index 
^=(2.8+1.3) between a maximum and a minimum scale of fluc- 
tuation of A max =(35±28)kpc and A m i n =(0.7±0.1)kpc, respec- 
tively. The corresponding magnetic field auto-correlation length 
is A B =5.2kpc. In addition, by including in the modeling the 
presence of X-ray cavities in coincidence with the radio galaxy 
lobes, we found a magnetic field strength of (Bo)=(l 1.7+9. 0)yuG 
at the cluster center. Farther out, the field decreases with ra- 
dius following the gas density to the power of ^=(0.9+0.5). To a 
comparison with other values reported in the literature, the radi- 
ally averaged magnetic field strength calculated over the central 
1 Mpc 3 is ~0. 19//G. The additional data and the numerical mod- 
eling of the intracluster magnetic field fluctuations allowed us to 
improve upon the previous estimate of 15//G (Eilek & Owen 
2002), with more stringent constraints not only on the magnetic 
field strength but also on its structure. 

From the literature the central magnetic field strength found 
in merger galaxy clusters is a few yt/G (e.g., 2.5//G in A2255 
by Govoni et al. 2006 and 4.7 yuG in Coma by Bonafede et al. 
2010). In cooling-core galaxy clusters the magnetic field central 
strength is still uncertain. Typically it is a few tens yuG (as in 
Hydra A, where values between 19 and 80//G have been indi- 
cated, e.g. Laing et al. 2008 and Kuchar & Ensslin 2011), even 
if lower values have also been found, as in the galaxy cluster 
A2634 (3 yuG, Vogt & Ensslin 2003). In A2199 we found a mag- 
netic field strength of about 10//G that seems to confirm val- 
ues of the magnetic field central strength in cooling-core galaxy 
cluster higher than in the merging cluster. 
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Fig. 9. Bayesian 3 -dimensional analysis of the RM structure 
function for the model with A max =35 kpc. Top panel: The 
dots represent the data (error bars are comparable to the size 
of the symbols). The shaded area represents the population 
of synthetic RM structure functions from the posterior dis- 
tribution. The dashed line corresponds to the most probable 
value for the model parameters (see text). Bottom panel: one- 
dimensional (histograms) and two-dimensional (colors and con- 
tours) marginalizations of the posterior for the model parame- 
ters. The contours are traced at 0.9, 0.75, and 0.5 of the peak 
value. 
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Fig. 10. The magnetic field power spectra (left panel) and radial profiles (right panel) for the models corresponding to the three 
values of A max =10, 35, and 164 kpc are shown as dotted, continuous, and dashed lines, respectively. The model with A max =35 kpc 
provides the best description of the observed RM structure function, see text. 
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Fig. 11. Qualitative comparison between the observed rotation measure image (left) and a synthetic 3-dimensional realization (right). 
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Appendix A: Maximum scale of fluctuation 

In this appendix we briefly discuss the determination of the outer scale of the magnetic field fluctuations. On the basis of the two- 
dimensional and three-dimensional simulations, we find that a maximum fluctuation scale of A max =35 kpc provides a very good 
description of the observed structure function. 

The observed rotation measure structure function is affected by a turnover in the large- separation regime (r ± » 2n/k m i n ). The 
drop comes from a lowering of the gas density at large distances from the cluster center, resulting in a systematic decrease of the 
power of the Faraday rotation on large scales. This effect cannot be modeled by the two-dimensional simulations (see top-left panel 
of Fig.0, while it is reproduced perfectly by the three-dimensional simulations (see top panel of Fig. [9]). However, since we know 
that the tapering imposed by the gas density distribution on the RM structure function at large separations may limit the possibility 
to determine the maximum fluctuation scale of the magnetic field, we may ask which is our actual sensitivity on A max . 

An accurate analysis would require considering A max as a free model parameter in the three-dimensional simulations as well, 
but the computational burden would be heavy. Indeed, we decided to explore only two more different values of A max , namely, a 
value of A max = 10 kpc, which is lower than the best-fit value, and a value of A max = 164 kpc, which is the maximum allowed by 
our computational grid of 1024 3 pixels. 

We performed the Bayesian analysis of the RM structure function with the same magnetic-field configuration as is described 
in §[5j except for the maximum scale of fluctuations. The results are presented in Fig. lA.ll for A max =10kpc and A max =164kpc (left 
and right panels, respectively). In the bottom panels of Fig. IA.ll we present the two-dimensional (colors and contours) and one- 
dimensional (histograms) marginalizations of the posterior. The magnetic field central strength and radial decrease are constrained 
to 

- <Bo>=(33.1±9.7)/iG, and t/=(1. 7+0.4) for A max =10kpc; 

- <£ >=(16.4±8.9)yuG, and ?7=(1.2±0.5) for A max =164kpc. 

In the top panels we show the observed structure function (dots) along with the synthetic structure functions from the posterior. The 
blue line is the best fit. We note that A max =10 kpc does not give a good description of the data because there is not enough power 
on large scales. 

On the other hand, a maximum scale of fluctuation of A max =164 kpc provides a better description of the observed structure 
function, even if it predicts slightly too much power for separations larger than r ± > 20 kpc. In this regime, the magnetic field model 
with A max =35 kpc is still better. 

Indeed, we may conclude that A max should be surely larger than 10 kpc. Although the drop in the gas density limits our sensitivity, 
we may infer that A max should be around 35 kpc, possibly lower than 164 kpc. 
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Fig. A.l. Bayesian 3 -dimensional analysis of the RM structure function for the model with cavities for A max =10kpc (left) and 
A ma x=164kpc (right). Top panels: The dots represent the data (error bars are comparable to the size of the symbols). The shaded 
area represents the population of synthetic RM structure functions from the posterior distribution. The dashed line corresponds to 
the most probable value for the model parameters (see text). Bottom panels: one-dimensional (histograms) and two-dimensional 
(colors and contours) marginalizations of the posterior for the model parameters. The contours are traced at 0.9, 0.75, and 0.5 of the 
peak value. 



